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We show that chiral kagome ice, an energetically stable manifold of Ising spins on kagome lattice, 
exhibits an anomalous quantum Hall effect when coupled to itinerant electrons. The highly degen- 
erate ice states provide a natural realization of power-law correlated random fluxes. Remarkably, 
the spectral gap of this system is robust despite the strong disorder experienced by the electrons. 
The presence of correlated disorder also stabilizes the a quantized Hall conductance over a wide 
range of filling fractions due to localization of electron eigenstates. We further show that emergent 
monopoles, i.e. topological defects in this ice-like manifold, bind a fluctuating electric dipole. 
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Integer quantum Hall effect is a cannonical example 
of the special role of topology in condensed matter. It 
is characterized by a topological invariant known as the 
first Chern number Besides the well-known two- 

dimensional electron gas in a magnetic field, quantum 
Hall effect can also emerge spontaneously from the inter- 
play of itinerant electrons and local magnetic moments 
in the absence of an external magnetic field @, Q ■ The 
origin of this phenomenon lies in the Berry phases im- 
parted to the electrons by magnetic textures exhibiting 
a nonzero scalar spin chirality Xijk = Sj • Sj x A 
periodic chiral order giving rise to a quantized Hall con- 
ductance can certainly be generated with long-range non- 
coplanar magnetic ordering: (Sj), (xijk) 7^ [4-1 No- 
tably, the scalar spin chirality can exist even without 
long-range spin order: (Sj) = and (xijk) 0, Hi- 
The discrete chiral order persists at finite temperatures 
even though the magnetic order is destroyed by thermal 
fluctuations in two dimensions. 

The stability of the Hall conductance in such metallic 
magnets is due to the robust local noncoplanar structure 
in a quasi-long-range magnetic order, as smooth distor- 
tions of the spin texture do not change the Berry flux 
pattern Q. However, it remains unclear whether a quan- 
tized Hall liquid can be stabilized in a state which does 
not exhibit magnetic ordering even at small length scales, 
giving rise to strongly disordered Berry fluxes. Here we 
demonstrate the existence of such exotic phase in a geo- 
metrically frustrated metallic magnet. More specifically, 
we show that a kagome-lattice model with electrons cou- 
pled to an extensively degenerate ice manifold exhibits 
a spontaneous quantum Hall effect. The quantum Hall 
ice phase considered in this work thus provides a proof of 
principle for a quantum Hall liquid coexisting with a clas- 
sical spin liquid. This work is partly motivated by recent 
experiments on a metallic spin-ice compound P^h^Cv, 
which shows an anomalous Hall effect in the absence of 
long-range magnetic order fl2- 14|. 

In a broader context, our work further illustrates 
the unusual electronic and transport properties resulting 
from correlated disorder 15, 3|. Despite great theoret- 



FIG. 1: (a) The kagome lattice and the Ising axes for 
local magnetic moments Si = en e;. (b) The projection of 
Si on the plane of the lattice for a random chiral kagome-ice 
configuration. Spins with a = +1 (a — —1) are shown in 
red (blue), (c) Two types of hopping amplitudes in a large- J 
spinless effective model, (d) The fluxes in the effective model 
through the hexagonal plaquettes for the ice configuration of 
panel (b) at the special canting angle 6* . 



ical interest in the effects of power-law correlated dis- 
order on conduction electrons, experimental realizations 
have remained elusive. Metallic magnets with geometri- 
cal frustration can provide such disorder quite naturally. 
For example, it was recently shown that the resistivity 
minimum observed in metallic spin ice can be explained 
by its peculiar spin correlations 17, HI- Here we show 
that a power-law correlated flux disorder is realized in 
the quantum Hall ice state. Although the quantized Hall 
conductance is expected to survive strong disorder, we 
also find that the spectral gap remains remarkably robust 
in the presence of this correlated random flux. In addi- 
tion to interesting bulk properties, topological defects in 
metallic spin systems also exhibit unusual phenomena 
such as charge fractionalization 11911 . Here we demon- 
strate that magnetic monopoles j20L |21| in the quantum 
Hall ice carry a fluctuating electric dipole, which could 
be detected in nuclear-quadrupole-resonance (NQR) ex- 
periments. 
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The specific model studied in this paper is a kagome- 
lattice system in which Ising-likc spins are subject to local 
constraints resembling the Bernal-Fowler ice rules 22 1. 
This so-called "kagome ice" [23| is an easy-axis ferromag- 
net with spins sitting on a two-dimensional network of 
corner-sharing triangles (Fig. [T|). The projections of the 
local easy axes on the kagome plane form a 120-degree 
ordering, while the axes are canted with respect to this 
plane by an angle 8. The direction of spins is specified 
by a set of Ising variables as S; = <Tjej. The magnetic 
charges (in natural units) for every up and down trian- 
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nearest-neighbor ferromagnetic exchange between spins 
Si can be recast into ^ Q Q 2 a I2i 



24l | . which penalizes 



triangles with magnetic charge ±3. It is thus energeti- 
cally favorable for each triangle to have magnetic charges 
±1. This implies the constraint that every triangle has 
either two incoming and one outgoing spins or vice versa. 

A subset of this kagome-ice manifold, known as the 
charge-ordered or chiral kagome ice 0, HBf , has a fur- 
ther constraint that spins in every up (down) triangle 
must be 2-in-l-out (l-in-2-out), i.e. all Qa = —1 and all 
Q v = +1 [Hi EH- Such configurations may be stabilized 
by two-body interactions of the form QaQv long-range 
dipolar interactions [13, HH , or alternatively in spin- ice 
pyrochlores subject to a magnetic field in the [111] di- 
rection (2(| [27j • Although there is no long-range order 
in the Ising variables cx^, the chiral kagome ice does have 
an overall magnetization in the out-of-plane direction, 
and our results do not directly explain the experiments 
on Pr2lr207, where the anomalous quantum Hall effect 
persists in the absence of net magnetization. 

We now introduce the itinerant electrons, which hop on 
the kagome lattice, and are coupled to the above Ising- 
like moments via the exchange coupling J. The electronic 
part of the Hamiltonian is then given by 



H = t E (4Sa+ H - 



a{ij) 



ia/3 



(1) 



where Ci a is a fermionic annihilation operator on site 
% and spin a. If the classical energy scales are much 
larger than the electron-mediated spin-spin interaction, 
coupling to itinerant electrons does not bring the sys- 
tem out of the ice manifold. It docs, however, select a 
particular configuration out of the chiral kagome ice man- 
ifold at zero-temperature. Our classical Monte-Carlo re- 
sults indicate that an ordered q = [Fig. |2{a)] state has 
the lowest energy, while the -\/3 x V3 [Fig. (2Kb)] is the 
highest-energy state. As we will see later, all chiral ice 
configurations give rise to a single-particle spectral gap 
A at 1/3 filling fraction, which is typically an order of 
magnitude larger than the difference, 5, between the en- 
ergy densities (per lattice site) of the q = and s/3 x \/3 
configurations. S is then the characteristic energy scale 
of electron-mediated spin-spin interactions. For example, 




FIG. 2: (a) The q = ice state is the ground state selected 
by electron-spin coupling at 1/3 filling, (b) The entropically 
favored y/3x v3 ice state has the highest energy when coupled 
to electrons. 



5 = 0.02i and A = 0.12t at 8 = 8* (a special angle to 
be discussed below). Therefore, the energy difference be- 
tween different chiral ice configurations is not resolved at 
temperatures 6 <C T <C min{ a } A({er}). For simplicity, 
we then work with an ensemble of Hamiltonians labeled 
by the Ising variables {a} at T = and 1/3 filling frac- 
tion. A quantized Hall conductance at T = thus im- 
plies a quantum anomalous Hall response at finite tem- 
peratures. We note in passing that even without QaQ\/ 
interactions, it may be possible to stabilize chiral kagome 
ice through coupling to itinerant electrons alone. 

To examine the intrinsic topological properties of the 
electrons in the ice manifold, we further simplify the 
problem by considering the strong coupling limit. We 
note that the spectral gap at 1/3 filling remains open 
even for finite J. In the |J| 3> \t\ limit, the electrons 
align themselves with the local moments and the effective 
hopping amplitude between two sites with local moments 
S 2 ; and Sj becomes t{xi\Xj)i where \xi) is the spinor 
cigenstatc of S, ■ <r a p ■ As shown in Fig. QJc) , there are 
two distinct hopping constants: t\ = cos f-e -1 ^" cos# and 
ti = sin^e 1 ^ + cos fe~*f sin£>, for opposite and same 
Ising spins on the bond, respectively. Note that these 
hopping amplitudes are written in a particular global 
gauge, while the electronic properties of the ice only de- 
pend on the gauge-invariant flux in each plaquette. The 
fluxes are equal in all up and down triangular plaquettes, 



4>i = arg(<i). 



(2) 



However, the fluxes in the hexagonal plaquettes of a 
generic chiral ice state are randomly distributed, giving 
rise to a strong off-diagonal disorder for the electrons. 
Nonetheless, the hexagonal fluxes are not uncorrelated. 
To see this, we note that two nearest-neighbor sites on a 
hexagon can not both have = —1, as each triangle in 
the chiral ice states contains one and only one such "mi- 
nority" Ising spins; see Fig.[IJb). Therefore, the number 
of bonds of type t\ is twice the number of sites with 
Cj = — 1. Thus 
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where X^ieO °* can take f° ur distinct values: 0, 2, 4, 
and 6. As we argue next, these fluxes are power-law 
correlated in the chiral kagome ice. 

The chiral ice manifold can be mapped to a dimer 
model on the honeycomb lattice, which connects the cen- 
ters of the corner-sharing triangles in the original kagome 
lattice 27 , [28 1 . By placing a dimer between the centers 
of two neighboring triangles connected by a minority spin 
<7j = — 1, each chiral ice state is mapped to a dimer cov- 
ering on the dual honeycomb lattice. The ensemble of all 
dimer coverings form a critical phase with dimer-dimcr 
correlations decaying asymptotically as 1/r 2 28| . As the 



dimers correspond to minority Ising spins, the linear re- 
lationship ([3]) implies that flux correlations obey 



(* O (r)*o(0)) - ($o(r)) 2 ~ 1/r 2 



(4) 



where (...) indicates an average over chiral kagome ice 
configurations. Note that the average flux in a hexagonal 
plaquette ($o( r )) = — ^4>\ ~ 202, is generically nonzero 
(the total flux through all triangles and hexagons van- 
ishes by construction). Numerically, we find that the 
flux distribution has a local y/3 x y/S pattern with a 1 /r 2 
decaying envelope. 

The magnitude of the hopping amplitudes can, in gen- 
eral, take two different values |ti| and \t2\- At a special 



canting angle 9 = 9* = ^ arccos (i), which we mostly fo- 



cus on in the present paper, we have |ii| = \t2\ = -4^ = t. 
For this special 9, the phases of the hopping amplitudes 
are given by 4>\ = — f and <p2 = fg, which according to 
Eqs. ([2]) and ([3]) leads to a flux — f in every triangular 
plaquette, and four different fluxes in the hexagonal ones, 
as shown in Fig. [TJd). 

In an ordered q = configuration, where the flux in 
all hexagons is equal to —2$ a, the tight-binding Hamil- 
tonian is known to exhibit a spectral gap and an integer 
quantum Hall effect at 1/3 and 2/3 filling fractions [4J. 
The topological origin of the q = ice can be understood 
by considering its band structure in the limits of 9 = 
and 9 = ^. The fluxes vanish in all plaquettes at 9 = 
but the tight-binding spectrum has pairs of Dirac points 
at 1/3 and 2/3 filling fractions. Non-zero fluxes resulting 
from canted spins gap out the Dirac points and lead to 
nontrivial band Chern numbers [1,13]. In the = f limit, 
we obtain an array of one-dimensional chains. Moving 
away from = 5 results in coupling these noninteract- 
ing Luttinger liquids (in the presence of time-reversal- 
symmetry- breaking fluxes), which also leads to quantum 
Hall effect [H,^. 

For random chiral ice states, the electrons experience 
flux disorder according to Eqs. (J3|) and (01 . Generically, 
flux disorder should not differ from electrostatic potential 
disorder if there is a net flux through the system [3l| . In 
case of the chiral kagome ice, the average flux vanishes, 
but since the system is characterized by a global quantum 
Hall response, one can expect similar behavior to quan- 
tum Hall systems, which emerge in the presence of a net 
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chiral kagome ice at 9 = 9* 




FIG. 3: (a) The density of states p(E) for chiral kagome ice at 
the special canting angle 9*. (b) The disorder-averaged Hall 
conductance as a function of the Fermi energy, (c) The ratio 
of density of extended states to the total density of states. 



magnetic field (note that the Hall conductance vanishes 
in a random magnetic field with zero average). 

Disorder often comes from the presence of impurities, 
which generically give an electrostatic potential with un- 
correlated fluctuations. Despite great theoretical inter- 
est, systems with power- law correlated disorder rarely 
occur in nature and need to be engineered [3J 



In 

our case, however, the relationship of the electronic Berry 
phases with a critical dimer model leads to a natural 
realization of power-law disorder. Generically, disorder 
closes the spectral gap but the quantum Hall effect nev- 
ertheless remains robust at T = 0. This can be under- 
stood in terms of the presence of a critical point cor- 
responding to a localization-delocalization transition be- 
tween two insulating states with ext ended states appear- 
ing at a single critical energy E c 



34 3 



For power- 
law correlated disorder with an exponent 2 as in Eq. (j4} , 
we expect this critical point to have a critical exponent 
v = 2.23 [32, 33 1, which characterizes the divergence of 
the localization length as \E — E c \~ v . As the transverse 
conductance stems only from extended states at a crit- 
ical energy E c (in the thermodynamic limit), filling the 
spectral gap with localized states does not change the 
quantum Hall response. 

While the quantized Hall conductance is expected to 
be robust, remarkably we find that the spectral gap at 
1/3 filling persists even in the presence of strong flux 
disorder. Fig. Eta) shows the disorder-averaged density 
of states obtained using loop-update Monte Carlo simu- 
lations at the special canting angle 9*. A spectral gap 
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A ~ 0.2t can be clearly seen. Interestingly, we also ob- 
serve a peak in the middle of the spectral gap. This 
peak corresponds to localized states around hexagonal 
plaquettes with $o = — 7r/2. The gap remains open for 
other canting angles, and the localized midgap states dis- 
appear. We have checked that the spectral gap (and a 
quantized Hall conductance) remains robust as long as 
the moments are noncoplanar (i.e., 9 ^ 0,7r/2). 

We also explicitly computed the quantum Hall con- 
ductance g xy using the real-space version of the Kubo 
formula a xy = Y^ m <T mi' where ^ indicates summation 
over occupied levels m, and 



a- = T Im 



(m\v x \n) (n\v y \m) 
(E m - E n f 



(5) 



|n) is a single-particle eigenstate with energy E n , A is 
the area of the system, and V{ is the velocity operator in 
direction i = x,y. We found that a xy is indeed quantized 
at filling fractions 1/3 and 2/3 for chiral kagome ice. 

Fig. EJb) shows the disorder-averaged func- 
tion of Fermi energy Ep. The Hall conductance rises 
much more quickly with increasing Ep for chiral kagome 
ice than the clean q = ice state, and reaches its quan- 
tized value at a smaller Ep. Similarly to in traditional 
quantum Hall states, where disorder leads to plateaus 
of transverse conductance as a function of the filling 
fraction, disorder stabilizes a zero-temperature quantized 
conductance over a wide range of filling fractions in chi- 
ral kagome ice. In order to suppress finite-size effects 
in the numerical calculation of cr xy (Ep), we average <r™ 
[Eq. (|SJ|] over two boundary phases <pi,2- For each real- 
ization of the disorder, the angle-averaged (ct™ ) is an in- 
tegral multiple of e 2 /h [3Sl - 44lj ] . A state |m) with nonzero 
( a xy) 7^ carries Hall current and is necessarily an ex- 
tended state [39[. 

Fig. [3(c) shows three peaks in the ratio of extended to 
total density of states. While the two higher energy peaks 
decrease with increasing lattice size, the lowest-energy 
peak persists and roughly coincides with the abrupt rise 
of cfxy' in Fig.^b). Despite the difficulty of precise finite- 
size scaling due to the large computation time required 
for both the disorder and boundary-angle averages, we 
estimate a critical energy E c = 1.9 ± 0.1, that is close to 
the bottom of the band. This is in contrast to the tra- 
ditional quantum Hall liquids, in which extended states 
appear in the middle of the broadened Landau level (the 
gap generally vanishes as Landau levels overlap). The 
proximity of E c to the bottom of the band implies that 
the Hall conductivity rises up more quickly to its quan- 
tized value in chiral kagome ice. 

We now consider the analogue of magnetic monopoles 
in the chiral ice manifold. These emergent excitations 
are topological defects violating the local constraints, 
which in the case of chiral kagome ice require Q a = — 1 
and Q v = 1. Emergent magnetic monopoles in the 




FIG. 4: (a) Flipping a string of spins (shown in red) in the 
q — state creates two magnetic monopoles. The variation 
of the magnetic flux in these triangles in compensated by a 
change of flux in the neighboring hexagon (shown in color), 
(b) Flipping many other loop results in a chiral kagome ice 
manifold with two magnetic monopole defects and global flux 
disorder, (c) The electric dipole in the distribution of electric 
charge in a q = state with two magnetic monopole defects, 
(d) The average distribution of electric charge for the chiral ice 
manifold with two monopoles indicating a fluctuating dipole. 
All data are for a system with 24 x 24 x 3 sites at filling 
fraction 1/3. 



chiral ice states arc thus defects with Qa = 1 and 
= —1 and carry magnetic charges ±2, respectively, 
with respect to a background of staggered charge or- 
der 



42|. Here the higher-energy defects that violate 



the ice rules (monopoles with charge ±4) are neglected. 
Two monopoles can be created in a q = state by flip- 
ping a string of head-to-tail spins as in Fig. HIJa). This 
process creates excess flux <f> and —(f) in a triangle and 
a neighboring hexagon. Through Laughlin's argument, 
one expects the adiabatic insertion of such flux to local- 
ize charges q = and — q in the corresponding triangle 
and hexagon, leading to an electric dipole, as indeed con- 
firmed in our numerical calculation shown in Fig. EJc). 

Since the Dirac string connecting the two monopoles 
is tensionless, these emergent monopoles are deconfincd 
in chiral kagome ice. The random fluctuations of the 
Dirac string in a generic chiral ice state gives rise to fluc- 
tuating dipolcs. It is worth noting that for individual 
realizations of disorder, the charge profile is highly dis- 
ordered in the bulk as various fluxes are inserted in the 
hexagons. Upon averaging, however, the charge varia- 
tions self-average in the bulk, and fluctuating dipole mo- 
ments with C3 symmetry emerge in the vicinity of de- 
fects (with a net quadrupole moment upon averaging). 
The averaged charge distribution around monopoles is 
shown in Fig. EJd), which clearly shows a C3 rotational 
symmetry. It may be possible to induce an average elec- 
tric dipole by breaking the C3 symmetry (e.g., through 
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changing the hopping amplitude of the horizontal bonds 
using pressure). 

In summary, the magnetic exchange between conduc- 
tion electrons and a chiral spin ice of localized Ising mo- 
ments leads to anomalous quantum Hall effect over a wide 
range of band filling factors. The critical correlations of 
the spin ice subsystem produce an effective power-law 
correlated disorder that has non-trivial consequences on 
the spectrum and transport properties of the conduction 
electrons. While previous studies have focused on the 
longitudinal conductivity of "metallic-ice" 



17, m 



have shown in this work that chiral spin-ice can dramat- 
ically change the electronic state by inducing a robust 
Quantum Hall liquid ("Quantum Hall ice"). Moreover, 
the interplay of the electrons with magnetic monopole de- 
fects of the ice manifold background leads to fluctuating 
electric dipoles with C3-symmetry. 
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